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Abstract 

The expected operating scenarios of 5G pose a great challenge to orthogonal frequency division multiplexing (OFDM) which 
has poor out of hand (OoB) spectral properties, stringent synchronization requirements, and large symbol duration. Generalized 
frequency division multiplexing (GFDM) which is the focus of this work, has heen suggested in the literature as one of the 
possible solutions to meet 5G requirements. In this work, the analytical performance evaluation of MMSE receiver for GFDM is 
presented. We also proposed precoding techniques to enhance the performance of GFDM. A simplified expression of SINK for 
MMSE receiver of GEDM is derived using special properties related to the modulation matrix of GEDM, which are described in 
this work. This SINK is used to evaluate the BER performance. Precoding schemes are proposed to reduce complexity of GEDM- 
MMSE receiver without compromising on the performance. Block Inverse Discrete Fourier Transform (BIDFT) and Discrete 
Fourier Transform (DFT) based precoding schemes are found to outperform GFDM-MMSE receiver due to frequency diversity 
gain while having complexity similar to zero-forcing receiver of GFDM. It is shown that both BIDFT and DFT-based precoding 
schemes reduce peak to average power ratio (PAPR) significantly. Computational complexity of different transmitters and receivers 
of precoded and uncoded GFDM is also presented. 

Index Terms 

Multi carrier modulation, BER, GEDM, PAPR, MMSE, Precoded GEDM. 


I. Introduction 

The need for wireless links, from proximity to a large distance communication, has been on the rise since its origin without 
ever showing any signs of distaste to growth. Towards this, several applications with diverse needs have driven the development 
of technical solutions. Amongst many, the domain of public wireless communication has provided one of the greatest benefits 
and drive to the development of society. The success is powered by standardization, which led to mass production and adoption 
of technology due to cost effectiveness. The development of IG - 4G HI, ||2l have fundamental drivers as the need for speed 
(higher data rate and higher spectral efficiency) and the requirement of lower latency. The next generation, 5G is also expected 
to be driven by similar demands but with added needs such as flexibility. Some of the important applications operating 5G 
0-© are flexible radio (cognitive radio), tactile internet, machine type communication (MTC), great service in crowd, super 
real time services, super reliable communication, and others. OFDM is de-facto transmission technology for broadband wireless 
access due to the robustness to multipath fading, high spectral efficiency and ease of implementation but has limited capability 
in meeting the needs of new requirements of 5G El, I©. MTC demands relaxed synchronization requirement which is a 
limitation in OFDM as it is sensitive to frequency errors. Small symbol duration is needed for low latency applications and 
due to the cyclic prefix (CP), OFDM will lead to low efficiency. High OoB radiation of OFDM limits, it’s use for opportunistic 
use of spectrum and dynamic spectrum allocation. 

To meet above requirements. Filter Bank Multi Carrier (FBMC) ll9l- lfT^ . Constant Envelope-OFDM (CE-OFDM) ifTSl . 
GEDM iflrl . Unified Eiltered Multi Carrier (UEMC) ifTSl have been proposed, where each has it’s own special features. 
EBMC, whose concept dates back to 1967 1(161 . is currently being considered as a candidate for 5G waveform due to it’s good 
frequency localization capabilities as well as low OoB radiation. However, EBMC signal suffer from inter symbol interference 
(ISI) as no cyclic prefix (CP) is used. The EBMC signal stretches in time and hence, limits spectral efficiency gains. Hence, 
EBMC is not suitable for low latency applications ini. UEMC (proposed for low latency applications) has limited suppression 
of OoB radiation. In this work, we focus on GEDM, which is another competitive waveform. GEDM has been shown to be 
flexible 0, ifTSl . 1191 in terms of using time-frequency resources. It has good OoB radiation properties ll20l . It is also quite 
resilient to synchronization requirement (HD and has good spectral efficiency as it uses circular pulse shaping which reduces 
cyclic prefix length in frequency selective fading channel(ESEC). 

Although GEDM has attractive features, yet there remains sufficient investigations to be done before it can be adopted into 
practical communication systems. One of the fundamental differences between GEDM and OEDM is that the former uses 
non-rectangular pulse shape on each subcarrier where it was rectangular for the latter. This causes inter-carrier interference 
(ICI) ED. Therefore, receiver design is an important issue in GEDM. It was earlier quite simple in the case of OEDM which 
led to its mass adoption. GEDM receiver in ESEC has been proposed as a two-stage receiver ED. Erequency domain zero 
forcing (ZE) equalization is used to mitigate the effect of the wireless channel in the first stage ifT^ . In second stage, linear as 
well as non-linear receivers have been proposed to mitigate self-interference in GEDM systems EHj IS). Linear receivers of 
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GFDM such as Matched filter (MF), Zero Forcing (ZF) and Minimum Mean Square Error (MMSE) are proposed in ||27I . It has 
been shown that ME is not able to mitigate self-interference. It has also been shown in ||22l, ll^ that MMSE performs better 
than ZE, especially when signal to interference and noise ratio(SINR) is low. Eor high SINR, MMSE performance matches 
with ZE. Analytical expression for bit error rate (BER) for ME and ZE receiver over additive white Gaussian noise channel 
(AWGN) is given in lIT^. EOll. Analytical expression of SINR is an essential part of performance analysis of a system as it 
leads to the calculation of error probability and capacity. BER computation for MMSE receiver is still not available in the 
literature. 

Since implementation of MMSE receiver requires very high amount of complexity, Successive Interference Cancellation 
(SIC) receivers to mitigate self-interference in AWGN channel were proposed in ll24l . Double SIC (D-SIC) receiver introduced 
in ll24l cancels interference from two adjacent subcarriers of all symbols. D-SIC can cancel out self-interference entirely in 
AWGN channel but requires few iterations thus induces processing delay at the receiver. Eor usability of GEDM, low complexity 
signal processing techniques remains to be investigated. 

Peak to Average Power Ratio (PAPR) is an important issue in the multicarrier communication.We expect GEDM to have 
a higher PAPR than GEDM as the use of identical non-rectangular pulse shapes increase PAPR in multicarrier system ll25l . 
Authors in flEj have compared PAPR of GEDM with GEDM for unequal number of subcarriers. As far as we know, PAPR 
comparison of GEDM with GEDM for equal number of subcarriers has not been investigated. Apart from that PAPR reduction 
schemes also need to be studied. 

A detailed exposition of the product of the modulation matrix with its hermitian reveals some interesting properties which 
helps in performance analysis and developing precoding techniques. The details are described in this paper. A simplified 
expression for SINR of MMSE receiver is developed using the above mentioned special properties. It is also shown that 
interference plus noise values can be approximated as a Gaussian random variable. Analytical BER is computed in AWGN 
and ESEC using derived SINR expressions. 

A detailed complexity analysis of different schemes is presented in this paper which shows that D-SIC is not quite simple 
to implement in comparison to ZE receiver. To reduce the self-interference at the receiver, precoding techniques for GEDM 
are investigated in this work. A generalized framework for precoding based GEDM is developed. Based on the properties 
related to the modulation matrix as mentioned above. Block Inverse Discrete Eourier Transform (BIDET) based precoding is 
proposed. Performance of BIDET, Discrete Eourier Transform (DET) and Singular Value Decomposition (S VD) based precoding 
techniques are compared with uncoded GEDM. PAPR of precoded GEDM is compared with uncoded GEDM and GEDM as 
well. 

The rest of the paper is organized as follows. The system model is developed in Section [III Analytical BER for MMSE 
receiver is presented in Section [Bll The precoding schemes proposed in this work are presented in Section |IVl Results related 
to the proposed works are given in Section |V] Section |Vl] has the conclusion. 

In this work vectors are represented by bold small letters, matrices are represented by bold capital letters and scalers are 
represented as normal small letters. The operation a(n modulus 1) is written as a{n)i, Ijv represents identity matrix with 
order N and represents convolution operation. E is expectation operation and j = y/—l. 

II. System Model 

GEDM is a multi carrier modulation technique with some similarity to GEDM. We begin by considering a block of QAM 
modulated symbols d = [dp di • • • We assume that data symbols are independent and identical i.e. E[didi] — 

, V( and E[did*] = 0 when I ^ q. Let the total bandwidth B be divided into N number of sub-carriers where symbol duration 
T = ^ second and ^ Hz is the subcarrier bandwidth. In case of GEDM this leads to orthogonal subcarriers. Let the symbol 
duration, T, be one time slot. GEDM is a block based transmission scheme and we consider a block to have M such time 
slots. Hence, in one block there are N subcarriers x M timeslots = NM QAM symbols. The flow of operations, as described 
below, can be understood in the light of Eig. [T] 

A. Transmitter 

Since we consider precoding, let P be a precoding matrix of size MN x MN, which is defined in Section |IV] The data 
vector d be multiplied with precoding matrix P and we obtain precoded data vector d = Pd. The MN x 1 precoded data 
vector d = [dp^p • • • dk,m ■ ■ ■ dN-i,M-i\^, where fc = 0 • • • — 1 denote subcarrier index and m = 0 • • • M — 1 indicates 

time slot index. Conventional GEDM system im can be seen as a special case of precoded GEDM system when, P = Iatm- 
The precoded data vector, d, is modulated using GEDM Modulator. Precoded data dk,m is first upsampled by N, which is 
represented as. 


= dk,mS{n - mN), n = 0,1, • • • , MN - 1. 


(1) 
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Fig. 1: Transmitter and receiver architecture for Precoding based GFDM system 


Now, this upsampled data is pulse shaped. Impulse response of pulse shaping filter is represented by g{n). Its length is MN. 
The upsampled precoded data is circularly convoluted with such pulse shaping filter g{n) and can be written as, 

MN-l 

= X! dk,mg{n-r)MN5{r-mN) = dk,mg{n - mN)MN- (2) 

r—0 


Now the filtered data is up-converted to fc**' subcarrier frequency and is given as. 


j2-Kkn 


^ — dk ,m{p) I 


(3) 


where, ak,m{n) = g{n — mN)MN& " for = 0, 1, • • • MN — 1, is called the kernel of GFDM for subcarrier and 
time slot. Now, output samples of all subcarrier frequencies and time slots are added and the GFDM signal can be given as. 


M-lAf-l M-lAf-1 

= X! X! ^k,m{n) = X! X! ^k,mak,m{jl). (4) 

m=0 k—0 m—0 k—0 

If we let, I = mN + k for m = 0, 1 • - ■ M — 1 and A: = 0, 1 - • • — 1. Then, we may write Q as. 


MN-l 

xin) = y] a^irip, (5) 

i=o 

where, superscript N is identifier of the specific mapping between index I and time slot and subcarrier index tuple {k, m) 
which takes N as scaler multiplier in the mapping. If we collect all output samples of GFDM signal in a vector called 
X = [x(0) a::(l) • • • x{MN — 1)]"^ and all samples of af^(n) is a vector called = [af^(O) af^(l) • • • af {MN — 1)]. Then 
using 0, X can be written as, 

MtV-l M7V-1 

X = df = df X column vector = Ajvd^, (6) 

1=0 1=0 

where, Ajy = [a^ ■ ■ ■ a^jy_j] is modulation matrix , = [do di ■ ■ ■ di - ■ ■ d^M-if^ is precoded data vector. If we 

collect all samples of ak,m{n) is vector called ak^m = [afe.m(O) afe,m(0) • • • ak,m{MN — 1)]"'’, then Ajv can also be written 
as, 

freq. freq. freq. 

An = [ao,o ai^o ■ • ■ ajv-i,o | ao,i ai^i • ■ ■ ajv-i,i | | ao,M-i ai,M-i ■ ■ ■ ajv-i,M-i] 

l®^time slot 2^‘^time slot Aff^time slot 


■ 9(0) 

9(0) 

. 9(0) 

p(mn-n) 

5{MN-N) 

... 

9(1) 

9 ( 1)9 

. 9(1)9”-^ 

p(mn-n+i) ••• 


(M-2)N terms 

.S(MN-l) 

S(MN~i)p“-l 

. g(MN-l)p'''-‘' 

^(MN-N-l) ••• 

■■■ 
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j2TT 

where, p = . At this point, it is also interesting to look into the structure vectors which constitute columns of A^r. 

The first column of A at i.e. aop holds all coefficient of pulse shaping filter and other columns of A^r or other vectors in the 
set of afc m’s are time and frequency shifted version of ap o where frequency index k denotes ^ shift in frequency and time 
index m denotes m time slot or mN sample cyclic shift. Taking clue from description of column vectors a^ structure of 
An can be understood by 0. Columns having same time shift are put together. Columns which are N column index apart 
are time shifted version of each other. Columns having same time shift are arranged in increasing order of frequency shift. We 
will explore time domain and frequency domain behaviour for an example An with total subcarriers N = A and total time 
slots M = 5. Pulse shaping filter is taken to be root raised cosine (RRC) with roll of factor (ROF) of 0.9. Fig. |2] shows the 
absolute values of each column index u = 0, 1 • • • MN — 1 with sample index n = 0, 1 • • • , MN — 1. It can be observed 
that first N = 4 columns have 0 time shift and next N columns have unit time shift and so on. 


_ 0.5 
■H 0.4 
c 0.3 
0.2 

“ 0.1 


20 



sample index n 


column index, u 


Fig. 2: Time Domain view of columns of A for N = A, M = 5 and ROF = 0.9 
Alternatively, if we take I = kM + m in equation 0, modulation matrix Am can be represented as, 

time time time 

Am = [a^ a^-• • ai,o ai,i ■ • • ai,M-i I ■ • ■ I ajv- 1,0 aAr- 1,1 ■ ■ • ajv-i,M-i]. (8) 

'-V-' '-V-' '-V-' 

isifrequency 2“‘^frequency iV^^frequency 

Matrix An and Am can be related as An = C-^m, where ^ is a permutation matrix which permutes column of matrix 
applied. Cyclic prefix is added to GFDM modulated block to prevent inter-block interference in FSFC. CP of length Ncp is 
prepended to x. After adding CP, transmitted vector, Xcp, can be given as 

x,p = [x(MA -N,p + 1: MN) ; x] (9) 

In the rest of the paper, for equations which are valid for both Am and An, modulation matrix will be denoted as A. However, 
wherever required An and Am will be specified. 
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B. Receiver 

Let, h = [/ii, /i2, • ■ ■ be L length channel impulse response vector, where, hr, for 1 < r < L, represents the complex 
baseband channel coefficient of r*** path ll27l . which we assume is zero mean circular symmetric complex Gaussian (ZMCSC). 
We also assume that channel coefficients related to different paths are uncorrelated. We consider, Ncp > L. Received vector 
of length Ncp + NM + L — 1 is given by, 

ycp = h*Xcp + Ucp, ( 10 ) 

where Ucp is AWGN vector of length MN + Ncp + L — 1 with elemental variance The first Ncp samples and last L — 1 
samples of yep are removed at the receiver i.e. y = [ycpiNcp + 1 : Ncp + MN)]. Use of cyclic prefix converts linear channel 
convolution to circular channel convolution when Ncp > L m. The MN length received vector after removal of CP can be 
written as, 

y = HAd + iz, (11) 


where H is circulant convolution matrix of size MN x MN, which can be written as. 


hi 0 
h2 hi 


0 hp ■ ■ ■ h2 

0 0 • • • /l3 


H = 


hp hp-i 
0 hp 


0 

0 


0 0 hp . hi 


( 12 ) 


and V is WGN vector of length MN with elemental variance a^. Received vector y is distorted due to (i) Self interference 
as subcarriers are non-orthogonal ll2^ and (ii) Inter-carrier interference(ICI) due to FSFC. 


III. Bit Error Rate Computation For MMSE Receiver 

A. MMSE Receiver 

To equalize the channel and GEDM induced self interference, a joint MMSE equalizer IT^ is considered here. Equalized 
data can be given as, 

d = [Imjv— + (HA)H(HA)]-i(HA)Hy 

(13) 

= Bd -t- Cu, 

where, B = [Imat^ + (HA)H(HA)]-1(HA)^^(HA) and C = + (HA)H(HA)]-1(HA)H. Eirst term in above 

equation holds desired plus interference values and second term holds the post processing noise values. 


B. SINR Computation 

Suppose we want to detect symbol. Estimated symbol can be given as, 

MN-l MN-1 


dl = [B]l,ldi+ 2 ^ [B]l,rdr 

r—0,r^l 


-f 


r—0 


[C]l,r 


(14) 


where first term is desired term, second term is interference term and third term is post processed noise term. Using above 
equation, E[\di\'^] = + cr^| where first term is average signal power Psig,i, 

second term is average interference power Psig+inr,i and third term is average post processing noise power Pnpp,i- Using 
([I3]i, E[dd^] = a'^BB^ + a^CC^, where diagonal values of first matrix term holds average signal plus interference 
power Psig+inr,i, diagonal values of second matrix term holds average post processing noise power, BB^ = [[Imat-^ + 
(HA)H(HA)]-1(HA)H(HA)]2 =02 and CC^ = [Ima^ + (HA)H(HA)]-1(HA)H(HA)[W^ + (HA)H(HA^]-i . 
Using this, average signal power Psig,i and average interference powerPinr,!, average signal plus interference power Psig+inr,i 
and post processed noise power Pnpp,i for F*' symbol, can be given as. 


Psig,l — 7 Psig+Inr,l — *^(^2 | [BB Plnr,l — Psig+Inr,l Psig,l Ultd Pjgpp^l — ]/^/| 

SINR for Z**' symbol can be computed as, 

7/ = 


PSig,l 


Pinr.l T Pnpp,1 

Now we will compute SINR for ESEC and AWGN channel separately. 


(15) 


(16) 
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1) FSFC: Both B and CC^ involve (HA)^(HA). To compute (fT^ we explore the product. (HA)^HA is a Hermitian 

matrix and hence it can be diagonalized as, (HA)^HA = VAV^, where, V is a unitary matrix which holds eigenvec¬ 
tors of (HA)^HA in it’s columns and A = Ai, • • • Amn-i} is a diagonal matrix which holds eigenvalues of 

(HA)HHA. We can write, B = VAVH, where, A = (img{Ao, Ai, • • • Amn-i}, As = and CC^ = VAV^^, where, 

A = diagjAo, Ai, • • • Amn-i} and As = — 7 ^^- . Putting values of B and C into (fTsT l. we can get, 

MN-1 

Ps^g,l=OdA[^kY%^l\^, Pir,r,i = ad 2 \ ^ [Vand PMpp,i = a, 2 \[vkv%,i\k (17) 

r—0,r^l 

Through the reduction, complex expression involving matrix inverse as in (fT4l i. is brought to a simpler form (using eigenvalue 
decomposition of (HA)^HA) i.e. instead of computing B we can proceed directly with HA. 

2) AWGN: In case of AWGN, H = Imtv, hence we can write the following, 

B=[^l + A^^A]-iA^^A, C = [^l + A^AJ-^A^ and CC“ = [^I + AHA]-1A^^A[^I + A^^A]-^ (18) 

d d d d 

It can be seen that A^A is a major component of the analysis. It is important to study the A^A before we proceed further. 
We will see the properties of Aj^Ajv below. 


AJ^Atv 

Matrix, An, can be decomposed as 

Atv = [[PoEjMATxAf [Pm-iEJmatxtv] , (19) 

where, Op’s be a set of MN x MN matrices, where p = 0, • • • ,M — 1. Suppose the first matrix in the set Oq = diag{g^} , 
where g = [( 7 ( 0 ) g{l) ■ ■ ■ g{MN — 1)]"'’. Any matrix in the set can be written as a circularly shifted version of Qg along it’s 
diagonals i.e. matrix in the set can be written as Op = diag{circshift[g^, —pN]}, where circshift represents right circular 
shift operation. Op can be described as. 


Qp = 


g{—pN)MN 


g{—pN -\- MN — l)MAr. 


( 20 ) 


J MNxMN 


and, E = [Ww ••• M times ••• Wat] is a MN x N matrix, where Wat is N x N normalized inverse DFT matrix 
[A^AAr]r -,5 = 0^0qE, where r,q = 0, - ■ ■ ,M — 1 and is N x N matrix. [A^AAr]r,g can be written as. 


[AjvAAf]r,<j(a,/3) = < 


a = j3 


J2k=o ^ otherwise 


( 21 ) 


where w = + K)MAr))g((—(ZiV + k)mn)), r,q = 0 ■ ■ ■ M —1 and a, (3 = 0 ■ ■ ■ N —1. Let, s = 0 • • • A — 1 

MN-l 

[AlANl^^{{a-s)N,{P-s)N)= Y. =[A^AN]r,q{(^,f3) ( 22 ) 


K—0 


This proves that each block [A^AAr]r,g of A^Atv is circulant. Let, = 0, • • • , M — 1. 

MN-l 

= Y + qN + K)MN)g{{-qN + qN + k)mn)- 


(23) 


K =0 


As w is periodic with N and p(/3) is periodic with MN, [A]^AAr](^_^)j^^ (a,/3) = [A]^Aa^]^ _^(a,/3), hence, A^^Aat 

is block circulant matrix with circulant blocks(BCCB) with blocks of size A x A. In the same way it can be proved that 
A^Am is also BCCB with blocks of size M x M. 

Using the properties of BCCB matrix one gets that 1^ (i) Inverse of BCCB matrix is a BCCB matrix, (ii) addition of a 
diagonal matrix with equal elements to a BCCB matrix is a BCCB matrix, (iii) multiplication of two BCCB matrix is BCCB 
matrix. Hence, it can be easily proved that B and CC^ is also a BCCB in case of AWGN channel. 

For any BCCB matrix Il29l it can be shown that (i) diagonal elements are identical , (ii) all rows have equal power and (iii) 
all columns have equal power. Using this it can be then concluded that, Psig,i = Psigi Psig+inr,i = Psig+inr, Pinr,i = Pinr 
and Pnpp,i = PNpp, VL Therefore we can write, [B]i^i = ^ Using this 
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and < fT5] ), we can write, 


Psig = Psig+Inr = -^iroce{B^} and Pnpp = -^trace{CC^}. 


Now, using (fTsT i in above equation, 

2 MN-1 . 2 MN-1 . 

= TT^^I £ -- Ps^g+Inr = ^ ^ 


MN-1 


s |2 


^0 A. + g 

SINR can be computed as. 


7 = 


s=o As + 


and P/vpp = 




E 


=0 (As + 


sig 


P, 


Inr+Pffp 


(24) 


(25) 


(26) 


Hence, SINR can be computed using the eigenvalues of A^A. Through the above, because of the BCCB property, we can 
compute SINR easily than by using the direct form. Inversion of matrix ^1 +A^A needs complexity of 0{M^N^) whereas 
eigenvalue computation of A^A needs complexity of 0{NM log 2 N) lE^ . Hence, SINR computation using above method is 
much simpler than using direct matrix computation. 


C. BER Computation 

1) FSFC Channel: Fig. shows cumulative distribution plot of signal plus interference value for 4 QAM modulated system 
for 3 dB and 9 dB ^ values. CDF plot for both cases is compared with Gaussian CDF with measured mean and variance 
values. It is clear from the hgure that interference plus noise values closely follow Gaussian distribution. Therefore the BER 




(a) FSFC channel 


(b) AWGN channel 


Fig. 3: CDF plot of interference plus noise value in AWGN Channel and FSFC for 4 QAM (Pulse shaping filter is RRC with 
ROF =0.5) 


for r'* QAM symbol with M. modulation order can be computed as 

v7M/2-1 


by 


Pb{E\-ii) ~ 4 


y/M-l 

VMlog2{M) 


E 


Q{{2r + 1) 


11 


(M-l) 


(27) 


7 / is the post processing SINR for symbol at the receiver given in (fTbl l. Average probability of error can be found as 


^ MN — 1 „cso 

X E / Pb(E\li)fiiili)dli, (28) 

where, is probability density function of 7 ;. 

2) AWGN Channel: Fig. shows cumulative distribution plot of interference plus noise value (fTTll for 4 QAM modulated 
system for 3 dB and 9 dB ^ values. As in case of FSFC, interference plus noise values closely follow Gaussian distribution. 
BER for QAM symbol with AA modulation order can be computed as ll30l . 


Pb{E) ~ 4 


vTw - 1 
VMlog2{M) 


VM/2-1 


E 


Q{(2r + 1) 



(29) 

































where 7 is the post processing SINR given in ( l26l l. 


IV. Precoded GFDM System 

In this section, precoding schemes to enhance performance of GFDM are described. The first precoding scheme BIDFT is 
described in Section HV-AI The scheme is developed using special properties of (HA)^HA (detailed in Section HV-All ) and 
A^A (detailed in Section IIII-Blb . For this precoding scheme there can be two kinds of receiver processing namely (i) Joint 
Processing ; described in Section IIV-AII which equalizes channel and GFDM modulation matrix simultaneously, whereas (ii) 
Two Stage Processing ; described in Section ITV-A2I first stage equalizes for the channel and second stage equalizes for GFDM 
modulation matrix. Two types of precoding matrices are defined for BIDFT precoding namely (i) BIDFT-M ; where A is 
structured in blocks of M x M i.e. Am and (ii) BIDFT-N : where A is structured in blocks of N x N i.e. An- BIDFT-N 
precoded GFDM is processed using Joint Processing as well as Two-Stage Processing whereas BIDFT-M precoded GFDM is 
processed using Two-Stage Processing only. 

DFT-based precoding is described in Section IIV-BI SVD based precoding is described in Section IIV-CI For each precoding 
scheme. Precoder matrix P, corresponding receivers and post processing SNR is described in detail. Both BIDFT and DFT-based 
precoding does not require channel state information (CSI) at the transmitter to compute P whereas SVD based precoding needs 
CSI at the transmitter to compute P. Channel knowledge at the transmitter can be maintained via feedback from the receiver 
or through the reciprocity principle in a duplex system llJll . BER for precoded GFDM system is presented in Section HV-DI 
Computational complexity of GFDM and precoded GFDM system is given in Section IIV-EI 


A. Block IDFT Precoded GFDM 

Received signal in (fTTTi can be processed in two ways, (i) Joint Processing ; Channel and self interference are equalized 
simultaneously, and, (ii) Two stage processing: Channel and self-interference and equalized separately. 

1) Joint Processing: Suppose the received signal in (fTTT) passed through a matched filter. Equalized vector which can be 
given as, 

MF ^ (HAP)Hy 

(30) 

= P“(HA)HHAPd + (HA)“iz. 

Since, (HA)^HA is multiplied to desired data vector in above equation. Now we will explore some properties of (HA)^HA 


(HA)^HA.- Using, description of An given in section ITlI-B II (HAAr)^HAAr matrix can be given as, 

[Lo,o]nxN ■■■ [Lo.M-l]ArxAf 


(HA^)"HA^ = 


[LM-l,o]AfxAf [LM-l,M-l]ArxAf, 


(31) 


MNxMN 


where Jjun = is a Ax iV sub-matrix or block, where T = H-^H is a MNx MN matrix and u,v = 0 ■ ■ ■ M —1. 

Channel convolution matrix H is a circulant matrix and it can be shown that Y = H^H is also a circulant matrix using 
the properties of circulant matrices 1281 . (HAAr)^HA 7 v will be a block circulant matrix, iff, = L„ „, where 

s = 0 ■ ■ ■ M — 1. In the expression of matrix and E are independent of block indices u, v. Therefore it can be 

said that, G is a block circulant matrix, iff, where, = Q^'EQv is MN x MN matrix. Let, 

= {vr,q}MNxMN and using the definition of Qp, it can be shown that 


— 9i—uN+r)MN9(—vN+q)MN'^’f'N and (32) 

(r, g^—(^u^s)N+r) MN9{ — {v+s)N+q)MN'^r,q- (33) 

Substituting, r' = r — sN and q' = q — sN and since, Y is a circulant matrix and hence V(^N+sN,q'+sN)MN ~ 'Ur',q'- Then, 
substituting r'& q' with r & q. 


^(u-l-siM,(«+«):«■ (^> 9 ) 9{-uN+r)MN9(-vN+q)MN'^r,q ^ (u) m ,(v) m i’'"’9) 


(34) 


Hence (HA 7 v)^HA 7 v is a block circulant matrix with blocks of size Ax A. (It may be noted here that (HAm)^HAm will not 
be block circulant with blocks of size M x M). Since (HAAr)^HA 7 v is block circulant matrix with blocks of size A x A, it can 

^ Wa '' 


r 

be decomposed as given in Il29l . l|32l, as, (HAAr)^HAAr = FbjvDf,jvF^^, where, F^n = [W 


Iat w]^1n 




-, where,W at = and E^n = diaglE'l 


where, WJy = --, wiicic,u;a' — c ^ anu ^j^n — '^‘‘’^yx^bN 

diagonal matrix with blocks of size A x A, where is diagonal matrix of size A x A. 
Using this decomposition and taking A as An, matched filter output in (l30l) can be written as, 

yMF ^ pHF.^Dfc^FH^Pd + (HA^)Hz.. 


'n ■ ■ 

. W^-l] 

N JmTVxMTV’ 

^In- 

• • is block 


(35) 

























9 


BIDFT-N Precoding: if we choose, P 
then, 

yMF._ 


Fbjv (call it block inverse discrete Fourier transform -N (BIDFTN) precoding). 


Dfeivd + (HAAr)Hiy 


\E>°bN ] 


r do^ 1 

DU 


df 

D^-i 

L ^bN J 


.<-1. 


+ i>, 


(36) 


where I> is MF processed noise vector. In above equation being block diagonal matrix, adds only — 1 interfering 
symbols instead MN — 1 (in case of uncoded GFDM). This shows that precoding reduces the number of interfering symbols 
significantly. Zero forcing equalization is applied to reduce interference further. Multiplying in above equation we get, 

d + (37) 

where, = D^j^(HAAr)^i>' is post processing noise vector and superscript JP signifies that signal processing steps followed 
in this method are joint processing (channel and self interference are equalized jointly). Df,jv can be computed as. 


Dbjv = pH^jHA^j^HAArF 


bN- 


(38) 


Post processing SNR for T'* symbol can be obtained as, 

'y-’P 

ibidft,l 






(39) 


where denominator in above equation is post processing noise power for symbol. 

2) Two Stage Processing: In above method, channel and GFDM were equalized together. In this method, we will first 
equalize channel distortions, and then GFDM induced self-interference. As explained in section |I] H is a circulant matrix. 
Hence, H can be decomposed as, 

(40) 


H = Watm^W 


H 

' NMy 


where, Wjvm is normalized IDFT matrix of size MN x MN and d/ = diag{vo, vi, - ■ ■ vmn-i} is a diagonal matrix. Channel 
equalized vector can be obtained as, 

=APd + WNM^-^W^MU, (41) 


where first term is the transmitted signal which is free from channel distortions completely, second term is enhanced noise and 
subscript FDE is acronym for frequency domain equalization (as above described channel equalization is frequency domain 
equalization lH]). ^ can be obtained as '5' = WJ^^HWatm which equivalently obtained by taking NM point FFT of zero 
padded channel convolution vector h which is also the first column of H HI, ||2l . 

Now passing channel equalized data Yfde to matched filter receiver, we can get, 

YFDE-MF = {AP)^yFDE = P^AHaPcI + (AP)HWwm’®'~^WH (42) 


It has been proved in section lni-B2l that A^A is BCCB matrix with blocks of size either N x N or M x M, which depends 
on whether A = Aat or A = Am- Which are described next. 

BIDFT-N Precoding: When modulation matrix is defined as Aat, A^An can be decomposed as, 

A^Aat = FbAfDftAfF^, (43) 

where, T>bN = diag{'D^^, is MN x MN block diagonal matrix where is r**' diagonal block of size 

N X N. Choosing, P = F^Af (BIDFTN precoding) and using above decomposition, matched filter output in (l42l i can be written 
as, 

yFDS-MJ’ = DbAfd + (AFfeArj^WATM'®' (44) 


Now multiplying in above equation, 

O-FDE-MF-ZF 


d + u 


N 

FDE-MF-ZF^ 


(45) 


where, iz^de_mf_zf — D^j^(AFf,jv)^WjvM’®' i^ enhanced noise vector. Since Df,jv needs to be computed to 

obtain dpuE-MFzF’ ii i’® computed as, Df,Ar = FbAf A^AF^^r- Post processing SNR for symbol can be obtained as, 

lFDE-MF-ZF,l - TTTjv TTn jHT77> 

^[^FDE-ME-ZE\^FDE-ME-ZE) \l,l 


where denominator in above equation is enhanced noise power for Z**^ symbol. 
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BIDFT-M Precoding: Now, if modulation matrix is defined as Am, A^Am = F^M^bAfF^, where, 

Im w'^mIm • 


' bM 


= [w^ 


w 


w 


w-il 


, where, = 


r(M—l)x 

'^M ImI ,-2, 

where,'u;M = 


' M " M " M JmTVxMTV’ ' M ^ 

and Df,M = diag{T)‘lj^^ ^Im ''' block diagonal matrix with blocks of size M x M, where is diagonal 

matrix of size M x M. Using this decomposition, choosing P = F^m (block IDFT-M (BIDFTM) Precoding) and following 
same signal processing steps as in case of A^v, equalized data vector can be given as. 


aM _ ^ 

O-fdE-MFzF — Cl ' 


,M 


where, 


FDE-MF-ZF 


'FDE-MF-ZF^ (47) 

= is enhanced noise vector. Since 'DbM needs to be computed to 


obtain dpoE-MFzF^ b can be computed as. 


DhM — Ffjjvf A'^AF^^. 


Post processing SNR for symbol can be obtained as, 

.M 


Ifde-mf-zfj — 




^W¥dE-M F-ZFi^¥DE-M F-Zf)^\i,1 


(48) 


(49) 


where denominator in above equation is enhanced noise power for symbol. In summary, BIDFT precoding can be understood 
by Fig. ga] 


B. DFT Precoded GFDM 

DFT precoding has been used in OFDM systems HI, El. It has been shown that DFT precoding reduces PAPR significantly 
and is one of the optimum precoding matrices to reduce PAPR in OFDM systems ||33l. This motivates us to investigate DFT 
precoded GFDM for reducing PAPR in GFDM system. Suppose, Q is spreading factor of the system then DFT order/ size 
can be computed as, N^ft = We assume that Q divides N completely. Two subcarrier mapping schemes are considered 
in this work 1^ (i) Localized frequency division multiple access (LFDMA) and (ii) Interleaved frequency division multiple 
access (IFDMA). 

Precoding matrix P can be defined as P = P^Pc, where Pc is a block diagonal matrix with each block being a DFT spreading 
matrix and P^ is a permutation matrix which implements subcarrier mapping i.e. LFDMA or IFDMA. The precoding spreading 
matrix Pc can be written as, 

Ndft 


Pc = 


w 


NdF7 


(50) 


where, W Noft normalized DFT matrix of size N^ft x No ft- Permutation matrix, P^ for LFDMA is an identity matrix. 
DFT precoded GFDM system can be understood from Fig. gb] Precoded data vector is GFDM modulated using the modulation 
matrix A. Received signal can be equalized using conventional linear im, 1231 or non-linear Il24l equalizer. We will present 
here ZF receiver for DFT precoded GFDM. ZF equalized precoded data vector can be obtained as. 


d,/ = (HA)-V 


(51) 


= d + (HA)~V. 

Equalized data vector ddft-spread can be obtained as, 

^dft —spread — P ^ 

De—spreading Subcarrier De—mapping 

where, Vdft-zf = P^((HA))“^i>' is post processing noise vector. Post processing SNR for Z*^ symbol can be obtained as, 

(53) 


= P“dz/ = P" X P" xd^/= d + izd/i-z/, 


(52) 


ldft-ZF,l 




where denominator in above equation is enhanced noise power for Z**^ symbol. 


C. SVD Precoded GFDM 

The product of the channel matrix and the modulation matrix (HA) can be decomposed as, 

HA = USV^, 


(54) 
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(b) DFT Precoded GFDM 


Fig. 4; Block Diagram of BIDFT and DFT precoded GFDM System 


where U and V are unitary matrix and S is diagonal singular-value matrix i.e. S = diag{sQ, si ■ ■ ■ Sr, ■ ■ ■ smn-i}, where Sr 
is singular value. Then, y in (fTTT l. can be written as, 

y = USV^Pd -F u. (55) 


At transmitter, by choosing P = V (assuming ideal feedback channel) and multiplying both sides with U^, the estimated 
symbol vector can be written as, 

= Sd -F U^iy. (56) 

The estimated symbol is then given by, 

MN-1 

dr‘^ = sidi+ ^ 

q=0 


(57) 
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From the above, it can be seen that by using SVD based precoding, interference is completely removed by orthogonalization 
of HA without the need for matrix inversion, which is required in zero forcing (ZF) and minimum mean square error (MMSE) 
receiver. SINR for symbol, can be computed as. 


7 r" = ^is/i 


(58) 


D. BER Performance of Precoding Techniques 

Expression for estimated data symbols, given in( l37l45l47l52l57l i, are summations of desired data symbols and enhanced 
noise vector. The enhanced noise vector is weighted sum of complex Gaussian random variable for a given channel realization. 
Hence, enhanced noise is also complex Gaussian random vector for a given channel realization. BER for QAM symbol with 
modulation order over ESEC can be obtained as. 


Pb{E\-fi)^A 


s/M-l 

VMlog2{M) 


VMI2-1 


E 


Q{(2r + 1) 


7i 


(M - 1) V 


where, 7 ; is post processing SNR computed in (I39I46I49I53I581 I. Average probability of error can be found as 

AIN — 1 poQ 


Pb{E) = 


1 


MN 


E 

i=0 


-P&(£’|7i)/7i Mdli 


(59) 


(60) 


where is probability distribution function of SINR for symbol. 


E. Computational Complexity 

Precoding schemes for GEDM are proposed in previous subsections. It is utmost important to compute and compare the 
complexity of these systems as it is directly proportional to the cost of the system. Out of different mathematical operations 
number of complex multiplication is a significant contributor to computational complexity m- Computation complexity is 
computed in terms of number of complex multiplication required to implement transmitter and receiver of precoded and 
uncoded GEDM system. It has been assumed that modulation matrix A is known at the receiver, hence any matrix that is 
derived from A is also known to receiver, such as, A^, A~^, etc. All known receivers for uncoded GEDM are considered for 
complexity computation of uncoded GEDM and DET precoded GEDM. Minimum number of complex multiplication required 
to perform various matrix and vector operations is computed . Complexity of SVD precoded GEDM receiver is computed for 
two cases namely (i)SVD of channel is (useful when channel is static for multiple transmit instances) and (ii) SVD of channel 
is unknown. Complexity computation can be found following Table H] 

V. Results 

In this section, the results related to works descri bed in earlier sections are presented. Analytical evaluation of B ER of 
MMSE receiver with GEDM is given in Section IV-AI BER evaluation of precoded GEDM system is provided in Section IV-BI 
Complexity of different transmitters and receivers of GEDM and precoded GE DM is given in Section IV-CI Einally, PAPR 
of proposed precoding schemes is compared with GEDM and GEDM in Section IV-DI GEDM system with parameters given 
Table [n] is considered here. It is assumed that the subcarrier bandwidth is larger than the coherence bandwidth of the channel 
for ESEC. SNR loss due to CP is also considered for ESEC. 


Number of Subcarriers N 

128 

Number of Timeslots M 

5 

Mapping 

16 QAM 

Pulse shape 

RRC with ROE =0.1 or 0.5 or 0.9 

CP length Ncp 

16 

Channel 

AWGN and ESEC 

Channel Length Nch 

16 

Power delay prohle 

[10~^]^, where a = 0,1 ■ ■ ■ Nch ~ 1 

Sub-carrier Bandwidth 

3.9 KHz 

RMS delay Spread 

4.3 fi sec 

Coherence Bandwidth 

4.7 KHz 


TABLE II: Simulation Parameters of GEDM receivers 


A. BER evaluation of GEDM with MMSE receiver 

BER vs ^ for 16-QAM in AWGN and ESEC for MMSE receiver is presented in Eig. |5] Legends ‘GEDM without Self- 

° 2 

Interference’ are obtained by using 7 ; = ^ in (l28T l for flat fading and in (|29] | for AWGN channel. This is a reference BER 
result for GEDM with zero self-interference and 16 QAM modulation in flat fading (l28T l and AWGN ( |29] |. The legends marked 
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Technique 

Operations 

Number of Complex 
Multiplications 

GFDM Tx 

one vector matrix multiplication for y = Ad 

(MNY 

GFDM Rx ZF- 
MF 

1) Frequency Domain equation (ZF) 

• one FFT operation for y' = Wy. 

• one diagonal complex valued matrix inversion and one diagonal matrix and 
vector multiplication for y" = A~^y'. 

• one FFT operation for ypDE = Wy'' 

2) Matched Filter for AWGN : one Matrix Vector multiplication for dzF = 
^^YFDE 

'^xog^MW) :r 

2MN + (MNf 

GFDM Rx ZF- 
ZF 

1) Frequency Domain equation (ZF): same operations as in ZF-MF. 

2) ZF for AWGN : one Matrix Vector multiplication for fde- 

Considering to be known and precomputed at the receiver. 

\og2{MN) ^ 

2MN + (MNf 

GFDM Rx ZF- 
MMSE 

1) Frequency Domain equation (ZF): same operations as in ZF-MF. 

2) MMSE for AWGN 

• MN complex value addition for computing C = ^ + A^A. 

• one matrix inversion for computing C~^. 

• one matrix vector multiplication for ytemp = A^y. 

• one matrix vector multiplication for dMMSE = C~^ytemp- 

(mnY 2{MNY + 

3 

GFDM Rx ZF- 
SIC 

1) Frequency Domain equation (ZF): same operations as in ZF-MF. 

2) DSIC for AWGN 

• MF operation to compute yMF MN(\/modorder — 1) comparators 

for detection of matched filter dieted ■ 

• 2M complex multiplication and 2M — 1 complex addition are needed for 
each sub-carrier and iteration index to compute yk = Adfc where k is 
sub-carrier index. 

• MN complex subtraction are needed for each sub-carrier and iteration 
index to compute yinterfree,k = YMF - Yk 

Xog^iMN) 

2MN + 2{MNYJ 

SVD precoded 
GFDM Tx 

two vector matrix multiplication for computing AVd 

2{MN)'^ 

SVD precoded 
GFDM Rx 

(Known SVD) 

one vector matrix multiplication for computing y 

{MNY 

SVD precoded 
GFDM Rx (Un- 
Known SVD) 

1) SVD computation 

2 ) one vector matrix multiplication for computing U^y 

(MTV)^ +26(M7V)^ 

BIDFT precoded 
GFDM Tx 

one vector matrix multiplication for computing AF^d as AF^ can be precomputed at 
transmitter 

{MNY 

BIDFT precoded 
GFDM Rx (Joint 
Processing) 

1 ) one vector matrix multiplication to compute HAF^ 

2) Computation of block diagonal matrix D 5 . 

3) Inversion which can be computed by inverting M square matrices of order 

N. 

{MNY \ogYN) =r 

2{MNY +mn'^ 

BIDFT precoded 
GFDM Rx (ZF- 
ZF) 

1) Frequency Domain Equalization : Same as in ZF-MF. 

2) AWGN Processing : 

• one vector matrix multiplication for computing ymf = {■^^b)^Y- 

• Computation of block diagonal matrix E)^^. 

• N times square matrix inversion of order M and block diagonal matrix 
and vector multiplication to compute d^oFT = '^^^Ymf 

\og2{MN) ^ 

2MN + NM^ + 

{MNY + 

NM^ or MN^ 

BIDFTM BIDFTN 

DFT Precoded 
GFDM Tx 

(additional over 
GFDM Tx) 

1) additionally, MQ times Noft point FFT and Sub-carrier Mapping 

2) same operation as in GFDM Tx 

log 2 NoFT 

DFT Precoded 
GFDM Rx 

(additional over 
GFDM Rx) 

1) same operation as in GFDM Rx. 

2) additionally,MQ times N^ft point FFT and Sub-carrier Mapping 

3) same operation as in GFDM Tx 

log 2 NofT 


TABLE I: Number of Complex Multiplication of different techniques in GFDM 
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‘Analytical’ are obtained by first computing 7 ; using ( fTbl ) for FSFC and using ( |26] ) for AWGN and then using this 7 ; in ( l28b 
and ( |29] | to obtain average BER for FSFC and AWGN respectively. The legends marked ‘Simu’ are obtained using Monte 
Carlo simulations using parameters given in Table |II] It is seen that BER from simulation matches quite well with with the 
analytical for both ESEC and AWGN. The gap between GEDM curve (analytical and simulation) and GEDM with no self 
interference curve is due to self interference encountered in GEDM ET\ . Erom the above discussion, it can be concluded that 
the expressions developed in this work are useful in estimating the theoretical BER for MMSE based GEDM receiver in ESEC 
and AWGN channel. 



0 2 4 6 8 10 12 14 16 18 

E,/NoindB 


Eig. 5; BER vs ^ for MMSE GEDM Receiver over AWGN and Erequency Selective Channel for 16 QAM (RRC pulse 
shaping filter is used with ROE=0.5) 


B. BER Evaluation of Precoded Techniques 

BER of precoded GEDM system is evaluated via Monte-Carlo simulations. Simulation parameters are given in Table lU One 
thousand channel realizations are used to obtain The BER results. The legends GEDM-ZE, GEDM-MMSE and GEDM-DSIC 
are used to represent the performance of the corresponding receivers for GEDM. The legends SVD-Prec, BIDETN, BIDETM, 
BIDET-JP, LEDMA-ZE and lEDMA-ZE are used to indicate the result of SVD based precoding, block IDETN precoding 
with two stage processing, block IDETM with two stage processing, block IDETN with joint processing, DET precoding with 
lEDMA and DET precoding with LEDMA respectively. Legend OEDM-CP is used to indicate theoretical BER performance 
of OEDM-CP. Pig. | 6 ] shows BER vs ^ for ROE = 0.1 in AWGN channel and ESEC. The following observations can be 
made. Under AWGN, all schemes have similar performance. It may be noted that in AWGN, there is no SNR loss as CP 
is not required. The performance of GEDM is only slightly worse than OPDM. In the case of ESEC, it is seen that SVD, 
BIDPT-N based precoding has performance similar to OPDM, as ICI is low because of small ROE. The better performance of 
DPT and BIDPT-M precoded GEDM over OPDM can be attributed to frequency diversity gain, which can be understood from 
Pigure |7] In this figure, it can be seen that one symbol in BIDPT-M has larger frequency spread compared to IPDMA which 
is larger than LEDMA. Accordingly, BIDETM precoded GEDM is performing better than IPDMA precoded GEDM which is 
performing better than LEDMA precoded GEDM. 

Next we look at Pigure 0 which shows similar performance analysis for ROE 0.9, which indicates higher ICI. In case of 
AWGN, degradation can be observed as compared to Pigure | 6 ] We see a similar degradation for all precoding schemes in 
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Fig. 6: BER vs ^ for Precoded GFDM Receiver over Frequency Selective and AWGN Channel using 16-QAM with 
iV=128,M=5, ROF=Q.\ (RRC) 


FSFC. This is expected for ROF 0.9 as significant overlapping of pulses in frequency. 

It is clear that the performance of precoding schemes is sensitive to ROF. In this work we have only considered RRC pulse 
shape. It is shown in llT9l . Il20l that performance of uncoded GFDM is sensitive to pulse shape choice. Proposed precoded 
GFDM system will also be sensitive to pulse shape choice as the enhanced noise expressions in (36, 43, 46, 51 and 56) are 
function of modulation matrix A which is a function of pulse shape as given in (7). 

It can be concluded that precoding schemes are giving good performance over FSFC. BIDFTM precoded GFDM performs 
best among all precoding schemes. 

C. Complexity Computation 

Table |IV] compares the complexity of different transmitters and Table |V] compares the complexity of different receivers for 
two different application scenarios given in ifT^ . Table Hill shows values of N and M considered for different application 
scenarios. Complexity is computed in terms of number of complex multiplications. Complexities of precoded GFDM systems 
and uncoded GFDM systems are compared with OFDM system. Spreading factor, Q = 4, is considered for DFT precoded 
GFDM. GFDM receiver processing for DFT precoded GFDM is considered to be GFDM-ZF. Number of iterations for DSIC 
receiver is considered to be 4 ll24l . 



Tactile Intemet(TI) 

Wireless RAN (WRAN) 

[N, M) 

(128, 5) 

(16, 127) 


TABLE III: Values of N and M for different application scenarios mi 

It is observed that the complexity of uncoded GFDM transmitter is about 100 times higher than OFDM transmitter in 
case of TI and, around 50 times greater than OFDM transmitter in case of WRAN. BIDFT precoded and DFT precoded 
GFDM transmitter has the same order of complexity as uncoded GFDM transmitter. Complexity of SVD precoded GFDM is 
around two times higher than uncoded GFDM transmitter. Therefore, it can be concluded that there is no significant increase 
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Fig. 7: Frequency spread of one symbol for different transmission schemes 


in complexity for BIDFT and DFT precoded transmitter when compared to uncoded GFDM transmitter. However for SVD 
precoded transmitter complexity is doubled, which is also not a significant increment. 

GFDM-ZF receiver complexity is much higher than OFDM receiver for instance 50 times in TI scenario. SVD precoded 
GFDM receiver when SVD is known, two-stage BIDFT precoded receiver and DFT precoded receiver have around same 
complexity as GFDM-ZF receiver. GFDM-MMSE receiver has very high complexity for example it 2000 times higher than 
GFDM-ZF in TI scenario. BIDFT precoded receiver complexity, when SVD of channel is not known, is even higher, for 
instance; it is 100 times higher than the complexity of GFDM-MMSE. It is also important to note that GEDM-DSIC has 
higher complexity than GEDM-ZE, DET precoded GEDM and two stage BIDET precoded GEDM. 


CASE 

OFDM 

Uncoded GFDM 




TI 

4.4 X 10'" 

4.09 X 10^' 




WRAN 

8.1 X 10^ 

4.1 X 10'’ 





TABLE IV: Number of Complex Multiplications for different Transmitters 


CASE 

OFDM 

GFDM- 

ZF 

GFDM- 

MMSE 

GFDM- 

DSIC 

SVD- 

GFDM Rx 

(KNOWN 

SVD) 

SVD-GFDM 
Rx (UN¬ 

KNOWN 

SVD) 

BIDFT- 

GFDM 

JOINT PRO¬ 
CESSING 

BIDFT- 

GFDM 

TWO 

STAGE 

PROCESS¬ 

ING 

DFT-GFDM 

TI 

8.9 X 
10 ® 

4.19 X 
10 ® 

8.8 X 

107 

3.2 X 
10 ® 





4.2 X 10“ 

WRAN 

1.62 X 
10 ^ 

4.16 X 
10 ® 

2.8 X 
10 ® 

3.3 X 
10 ^ 

HljmH 

2.1 X 10^^ 

2.4 X 10' 

mm 

4.16 X 10® 


TABLE V: Number of Complex Multiplications for different Receivers 
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Fig. 8: BER vs ^ for Precoded GFDM Receiver over Frequency Selective and AWGN Channel using 16-QAM with 
iV=128,M=5, ROF=0.9 (RRC) 


D. PAPR of Precoding Techniques 

The impact of precoding on PAPR is presented in Fig. |9] Complementary cumulative distribution function (CCDF) of PAPR is 
computed using Monte-Carlo simulation. 10® transmitted blocks were generated, where each block has two precoded GFDM 
symbols. For each precoded GFDM symbol : N = 128, M = 5 and ROF for RRC pulse shaping filter = 0.5. For OFDM, 
iV = 128. Randomly generated data symbols are considered to be QPSK modulated. 

As expected, GFDM has worse performance than OFDM. Precoding has a positive effect on GFDM. We compare the PAPR 
value that is exceeded with probability less than 0.1% (PrjPAPR > PAPRo = 10“®}). SVD precoded GFDM and BIDFTN 
precoded GFDM reduces the PAPR by 0.3 dB, but it is still higher than OFDM. DFT precoded GFDM with LFDMA subcaiTier 
mapping reduces PAPR by 3.4 dB and is lower than OFDM. DFT precoded GFDM with IFDMA subcarrier mapping and 
BIDFTM precoded GFDM reduces PAPR by 9 dB. 


VI. Conclusion 

GFDM, a possible waveform for 5G with flexibility to access time-frequency radio resources, is investigated in this work. 
Two properties of GFDM system are proved (i) product of modulation matrix with its hermitian is a block circulant matrix 
with circulant blocks and (ii) product of circulant convolution channel matrix with the modulation matrix when multiplied 
with its hermitian is block circulant matrix. An expression of SINR for MMSE receiver for GFDM is developed using matrix 
representation for the signal model. The SINR expression is determined in terms of eigenvalues of the product of modulation 
matrix with its hermitian matrix for AWGN channel and in terms of eigenvalues and eigenvectors of product of circulant 
convolution channel matrix with the modulation matrix when multiplied with its hermitian for FSFC. It is also shown that the 
addition of interference and noise values are Gaussian distributed. Using this BER is found. It is found that the BER obtained 
from simulation and that from the expression developed here match quite well in AWGN channel as well as in FSFC. 

Three new precoding techniques which improve performance of GFDM system are proposed. The SVD based precoding for 
GFDM removes interference by orthogonalizing the received symbols. It does not require matrix inversion yet its performance 
is quite close to that of GFDM-MMSE result. The complexity of the receiver for such precoding is quite high when compared 
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Fig. 9: Complementary Cumulative Distribution function (empirical) of PAPR for GFDM for A^=128 and M=5, ROF =0.5 
(RRC) and Q= 4. 


to GFDM-ZF receiver when SVD of channel is not known. However, when SVD of channel is known complexity of the 
receiver is comparable to GFDM-ZF. Hence, SVD based precoding can be used in cases when channel is constant for multiple 
transmit instances. Performance of BIDFTN precoding is found to be better than GFDM-ZF and GFDM-DSIC. Two stage 
BIDFTN precoded GFDM has complexity similar to GFDM-ZF and lesser than GFDM-DSIC as well as gives lower PAPR. 
Hence, two-stage BIDFTN precoded GFDM should be preferred than GFDM-ZF. BIDFTM precoded as well as DFT precoded 
GFDM receivers which require complexity similar to GFDM-ZF performs much better than even GFDM-MMSE receiver under 
FSFC. Apart from this, BIDFTM and DFT precoded GFDM reduces PAPR significantly. Precoded GFDM system proposed 
in this work can give better BER performance than GEDM system with no increase in complexity with the added advantage 
of decreased PAPR. It can be concluded that BIDETM precoded GEDM should be preferred to other precoding schemes as 
it gives better BER performance than other precoded and uncoded receivers and decreases PAPR significantly without any 
increase in complexity. 
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